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STEADY-STATE HOT-ION PLASMA PRODUCED BY 
CROSSED ELECTRIC AND MAGNETIC FIELDS 

by Donald R. Sigman and John J. Reinmann 

Lewis Research Center 

SUMMARY 

Results of initial diagnostics on the Lewis Research Center hot-ion-plasma source 
(HIP-1) are reported. The mirror-contained plasma was heated by applying a radial 
electric field similar to that used in the Oak Ridge National Laboratory (ORNL) Burnout 
experiments. An electrostatic neutral particle analyzer was used to make a parametric 
study of ion energy distributions in both hydrogen and deuterium plasmas. Ion temper- 
atures as high as 2 keV were measured in plasmas with electron densities greater than 
12 -3 

10 cm . The performance of the device was found to be extremely sensitive to mag- 
netic field. There are indications that ion heating was reduced when the size of the 
Larmor orbit was larger than the region of strong radial electric field. 

INTRODUCTION 

This report describes the first experimental data obtained with the HIP-1 (hot-ion 
plasma) source at the NASA Lewis Research Center. In this facility, energy is coupled 
to the plasma by applying a radially inward dc electric field £ near the maximum field 
region of a magnetic mirror. The mutually perpendicular electric and magnetic fields 
cause the plasma to rotate (drift) in the azimuthal direction with a velocity c<?/B . 

(All symbols are defined in the appendix. ) According to the theory of Hirose and Alex- 
eff (ref. 1) an aximuthal current exists because of a drift velocity difference between 
ions and electrons. The combination of this current and a radial density gradient leads 
to the growth of electrostatic azimuthal waves at frequencies near the lower hybrid fre- 

/r / 2 / 2 \i 1/2 

quency, u> LH = / 1 + ( w p e / w ce J • Hirose and Alexeff then postulate that ions 

are heated in planes perpendicular to the applied magnetic field by the high-frequency 

fields of the wave. A similar experiment, the Burnout (ref. 2) experiment at the Oak 

Ridge National Laboratory (ORNL), has yielded ion temperatures as high as 3 keV and 

13 -3 

particle densities greater than 10 cm 



DESCRIPTION OF EXPERIMENT 


A schematic of the experimental apparatus is shown in figure 1. The maximum 
value of magnetic field at the midplane B Q was 1. 8 teslas and the mirror ratio was 
1. 58. The voltage available to be applied to the electrodes was 20 kilovolts. The gas 
feed was through the center of the cathode. The electrodes were identical to those used 
in the Burnout VI experiment at ORNL. 

There are several modes of operation of the facility, depending on the background 
gas pressure. At low pressures (below 5x10" 5 torr) a faint glow discharge exists, and 
the full 20 kilovolts of the power supply can be applied without arcing. However, less 
than 0. 1 ampere is drawn. At higher pressures (above 2xl0" 4 torr), very high currents 
exist (4 to 6 A); but the voltage between the electrodes drops severely as a result of a 

ballast resistor in series with the plasma. At intermediate pressures of approximately 

_ ^ 

1x10 torr, a high-voltage (8 to 10 kV), intermediate -current (1 A) mode exists. Pre- 
liminary ion temperature measurements indicated the highest values for this mode, and 
consequently all the data presented hereinafter are at p« lxlO -4 torr. 

The current -voltage characteristics for the intermediate-pressure mode are as fol- 
lows: As the voltage between electrodes is first raised from zero, very little current 
flows. Then at about 2 to 4 kilovolts the current starts to rise rapidly until a nearly 
saturated value is reached. Further raising of the voltage then leads to only slightly 
higher currents. Finally, a limitation is reached on the voltage when an arc occurs be- 
tween the electrodes. This arc overvoltage is found to be from 9 to 12 kilovolts, being 
higher for larger magnetic fields. The ballast resistor mentioned previously is used to 
limit the current when it first begins to rise rapidly. 

Conditioning was also found to be an important factor in achieving the hot-ion mode 
previously described. Whenever the system has been open to atmosphere or just sitting 
overnight, the high-current mode prevails when the voltage is applied. It was assumed 
that this is caused by a pressure increase due to outgassing from the carbon electrodes 
as they heat up. After several runs, however, the outgassing either ceases or reaches 
a constant level and repeatable operation in the hot-ion mode is possible. Hydrogen 
seems to condition more rapidly than deuterium. The highest currents also occur in 
hydrogen (up to 1. 2 A); whereas, the maximum achievable current in deuterium is about 
0. 8 ampere. Any attempt to raise the currents above these values by raising the gas 
flow rate always leads to the high-current mode. 


DIAGNOSTICS 

Most of the data herein come from two principal diagnostics: (1) the measurement 
of the energy spectrum of charge -exchange neutrals escaping the plasma and (2) the 
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measurement of the Doppler -broadened 486. 1-nanometer /3-line of hydrogen or deu- 
terium. 


Neutral Particle Energy Analyzer 

A schematic of the neutral particle energy analyzer (NPEA) and its arrangement 
with respect to the plasma is shown in figure 2. The NPEA was a duplicate of the one 
described in reference 3, but with a different particle detector. After emerging from 
the plasma, the collimated beam of neutrals was ionized in the nitrogen gas cell. The 
energy distribution of the beam of positive ions was obtained by sweeping the voltage on 
the 90° electrostatic deflection plates. The detector was a Bendix model 306 magnetic 
electron multiplier (MEM). The current output of the MEM was fed into a Keithley 
model 417 picoammeter. The output of the picoammeter went to the y-input of an x-y 
recorder. The NPEA sampled, from along a chord of plasma, a region approximately 
1 centimeter high (vertical) and 2. 5 centimeters wide (parallel to the magnetic field). 

The NPEA apparatus was mounted so that the plasma could be scanned vertically by 
pivoting the beam line about the point located at the test section flange so designated in 
figure 2. 

The ion energy distribution within the plasma was obtained from the energy distri- 
bution at the detector output in the following manner: The detector current is propor- 
tional to the product of several factors 

1(E) ~ ff cx (E)v An E k qc K EDP K D ^ 


where 


1(E) 

a cx< E > 

V 

An E 

Kgc ( e ) 

k edp^ e ^ 

K d (E) 


detector current originating from ions in plasma with energies between E 
and E + AE 

charge exchange cross section between ions and molecular hydrogen 

o 

speed of ions with energy E = 1/2 nijV 

number of plasma ions per unit volume with energies between E and E + AE 
neutral-to-ion conversion efficiency of nitrogen gas cell 
energy resolution of electrostatic deflector plates 

secondary emission coefficient for plasma ions incident on MEM tungsten 
cathode 
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The charge exchange cross sections cr were taken from reference 4. The gas 

cx 

cell conversion factors presented in reference 5 were used. The energy resolu- 

tion of the electrostatic deflector plates K EE p was proportional to E/AE. Least 
known was the secondary emission coefficient K D for plasma ions incident on the 
MEM's tungsten cathode. We assumed that K D is proportional to E. The variation 
of An E with energy depended on the energy distribution within the plasma. We assumed 
an isotropic distribution in velocity space which resulted in the following relation: 



In terms of energy, An E has the following proportionality: 

An E ~ ^E f(E)AE 

Substituting these energy dependences into equation (1) yielded the energy distribution 
function f(E) in the plasma 


f(E) ~ ^ (2) 

ct cx< E > e 3 k gC< E > 

Doppler Broadening 

To make the Doppler -broadening measurements we employed a simple two-slit col- 
limator to view along a chord of the plasma, the line of sight being horizontal and per- 
pendicular to the axial magnetic field. The light sampled came from a plasma region 
approximately 1 centimeter vertically by 7. 5 centimeters axially, and vertical scans 
could be made to look along various chords. The spectra were taken with a Jarrell-Ash 
0. 5-meter Ebert spectrometer. 

Other diagnostics used so far to take a cursory look at the plasma were 8-millimeter 
microwave interferometry for electron density, Langmuir probes for plasma potential, 
and relative intensity measurements of helium lines at 504. 8 and 471. 3 nanometers for 
electron temperature. 
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EXPERIMENTAL RESULTS 


Neutral Particle Measurements 

Figure 3 shows a comparison of the NPEA detector current as a function of energy 
for hydrogen and deuterium for a discharge with an applied voltage V Q of 8 kilovolts 
and a B Q of 1. 35 teslas. These are typical spectra for almost all values of V Q and 
B q for which the hot-ion mode exists. In deuterium, the peak was broader but the high- 
energy tail was always less relative to the peak than is the case for hydrogen. 

Figure 4 shows reduced energy spectra f(E) for the best achievable conditions of 
V Q = 10 kilovolts and B Q = 1. 8 teslas. In reducing these data it was assumed that the 
charge exchange was between molecules and ions. We see that the temperatures indi- 
cated by the slope of the curves at high energies are 2 keV for hydrogen and 1 keV for 
deuterium. The two nearly identical curves for deuterium indicate that the shape of the 
spectrum is roughly the same whether one is viewing directly through the center of the 
plasma or along a chord passing 1. 5 centimeters below the center. Figure 5 shows two 
of the curves from figure 4 replotted as a function of AE where A is 1 for hydrogen 
and 2 for deuterium. Equal values of AE correspond to equal values of the Larmor 
radius. 

As one scans the plasma and looks at successive chords both above and below the 
centerline of the plasma, there is evident an asymmetry in the relative amplitudes of 
the spectra. This is caused by the fact that in a rotating (drifting) plasma there are 
more particles of a given energy moving toward the analyzer which are below the center- 
line than there are above the centerline. A reversal of the magnetic field direction re- 
versed this asymmetry. It is unfortunate that this drift-related azimuthal variation in 
the velocity distribution makes it impossible to Abel-invert the data to determine the 
radial variation of the distribution function. 

Figure 6 shows reduced data (V Q = 8 kV, B Q = 1. 35 T) taken along a line of sight 
through the center of the plasma and also along a line of sight 1. 5 centimeters below the 
center. The qualitative aspects are the same for hydrogen and deuterium. The curves 
for each gas have been normalized at 1 keV. One possible interpretation of these data 
is that the temperature of the plasma is peaked near the center; hence, the tail of the 
distribution is enhanced when sighting through the plasma center. This was not found to 
be the case for the same two sighting positions at V Q = 10 kilovolts and B = 1. 8 teslas 
(fig. 4). 

It also may be that the highest energy particles were carbon (or some other impur- 
ities). A momentum analyzer, now being built, will tell us whether species other than 
H (or D) contribute to the measured currents. One must realize that a clear interpreta- 
tion of such data is difficult because (1) the effect of plasma drift on the distribution 
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changes with each new line of sight, and (2) the data presented are an integration along 
a chord of the plasma. 

Figure 7(a) shows hydrogen raw data spectra with B Q as a parameter and V Q 
fixed at 8 kilovolts, while figure 7(b) has V Q as a parameter and B Q fixed at 1. 8 tes- 
las. Part of these data are then replotted in figures 7(c) and (d) as the ratio 
I(E)/I(1. 0 keV) as a function of B Q and V , respectively. An increase in this ratio 
means there is a higher fraction of the ions at high energies, and we interpret this as 
ion heating. The data show an apparent increase in ion heating with increases in B Q 
and/or V . 

One can get an idea of the plasma diameter by cross plotting data at a given particle 
energy from spectra taken along various lines of sight. Figures 8(a) and (b) show such 
plots for both hydrogen and deuterium and for several energies. The mean plasma 
diameter was approximately 2 centimeters and was slightly larger for deuterium than 
for hydrogen. If it were possible to Abel-invert the data, the mean diameter would 
probably be slightly larger. The shift of the data peaks away from the centerline sight- 
ing position was again caused by the drift. For the higher energy particles, the Larmor 
orbits were large enough to cover a sizeable portion of the plasma. This tended to 
’’smear out" the drift asymmetry. 


Doppler Broadening 

Figure 9 shows typical H0 spectra (V Q = 8 kV, B Q = 1. 35 T) for sightings below, 
through, and above the centerline, respectively. The center of the broadened base, 
which comes from the fast charge exchange neutrals, shifts from below 486. 1 nano- 
meter to above as the line of sight is shifted upward. The maximum shift away from 
486. 1 nanometers is about 1. 2 nanometers, which corresponds to a plasma drift veloc- 
ity of 7. 4x10^ centimeters per second. Later we discuss measurements which indicate 
that the potential on centerline at the midplane was approximately 40 percent of that ap- 
plied to the cathodes. Using this value, we calculated an average electric field in the 
previously discussed annular region between the electrodes of 3 kilovolts per centime- 
ter. The c«?/B drift velocity was then computed to be 2x10^ centimeters per second, 
or about three times that measured by the shift. This is not necessarily inconsistent 

when one considers the fact that the lifetime of the excited states is of the order of 
-8 

10 second. Thus, the hot neutral atoms may travel a significant distance across the 
plasma before radiating, making the origin of the light quite different from the position 
of the original ion. This effect will tend to reduce the measured drift. 

When G(A) from figure 9 is plotted logarithmically as a function of (AX)^ (where 
(AX) ~ E), the broadened base part of these curves gives a straight line. The temper- 
ature related to the slope of this line is 0.35 keV for hydrogen. Raising the magnetic 
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field to 1. 5 teslas raised this value to 0. 40 keV. These non-Abel-inverted ion temper- 
ature values determined from the spectroscopic data probably represent a mean value 
for the whole plasma; whereas, the higher values determined from the tail of the NPEA 
data are either indicative of the peak temperatures or represent data for other species. 

Figure 10 shows the peak of the Hj3 curves plotted as a function of sighting position 
(V Q = 8 kV, B q = 1. 35 T). The plasma diameter indicated here is the same as that de- 
rived from the neutral particle data (~2 cm). 


Other Diagnostics 

The beam of the microwave interferometer used to make electron density measure- 
ments was larger than the plasma diameter. This is certainly far from an ideal condi- 
tion and would lead to a measured phase shift for the whole beam which is less than the 
phase shift of the portion of the beam which passed through the plasma. Thus, all we 
can say is that the mean density n g of the plasma is greater than 10^ cm - **. 

Electron temperature T g measurements using the helium line intensity ratio 
method were made only in hydrogen for V Q = 8 kilovolts and B = 1. 2 teslas. The 
measured value for T g was 12 eV. Since these are conditions of low power input to 
the plasma, it is probably safe to say that the electron temperature for the conditions 
discussed herein was greater than 12 eV. 

A probe actuator was used to inject a Langmuir probe to the center of the plasma 
and out again in less than 0. 1 second. At 4 and 6 kilovolts the probe measured potentials 
on axis of approximately 40 percent of the applied voltage. Under these conditions the 
probe became incandescent during the traverse. At 8 kilovolts the probe disintegrated. 
These measurements give further evidence that a large radial electric field exists in the 
midplane as well as at the mirror throats. 


CONCLUDING REMARKS 

Initial diagnostics on the HIP-1 hot-ion-plasma source revealed more than one mode 
of operation. Highest ion temperatures were observed in a mode where the total 
cathode-anode current was limited to approximately 1 ampere. Highest temperatures 
measured were about 2 keV in hydrogen and 1 keV in deuterium. Plasma heating was 
found to increase with applied voltage and/or dc magnetic field strength, being more 
sensitive to the magnetic field. Electron densities were greater than 10 13 cm -3 . Evi- 
dence of plasma rotation was seen through asymmetries in radial scans of the neutral 
particle energies and the Doppler -broadened H/3 line. Determination of the ion energy 
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distribution function was difficult because of the strong role played by ion drifts in de- 
termining the measured neutral particle energy spectrum. 


Lewis Research Center, 

National Aeronautics and Space Administration, 
Cleveland, Ohio, February 7, 1973, 
503-10. 
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APPENDIX - SYMBOLS 


A atomic number for ions 

B Q axial magnetic field at midplane 

c speed of light 

E ion or neutral particle energy 

& radial electric field 

f plasma energy or velocity distribution function 

G spectrometer current 

I neutral particle analyzer detector current 

Kj-j neutral particle analyzer detector efficiency 

Kedp electrostatic analyzer resolution 

K GC gas stripping cell conversion efficiency 

m i ion mass 

n g electron density 

p background gas pressure 

T g electron temperature 

V Q applied voltage between electrodes 

v ion or neutral velocity 

\ wavelength 

a charge exchange cross section for H (or D) ions on molecular 

u> ce electron cyclotron frequency 

o>Lh lower hybrid frequency 

w electron plasma frequency 

pe 

ojpj^ ion plasma frequency 



REFERENCES 


1. Hirose, A. ; and Alexeff, I. : Electrostatic Instabilities Driven by Currents Perpen- 

dicular to an External Magnetic Field. Nucl. Fusion, vol. 12, no. 3, May 1972, 
pp. 315-323. 

2. Anon. : Thermonuclear Division Annular Progress Report for Period Ending Decem- 

ber 31, 1970. Rep. ORNL-4688, Oak Ridge National Lab. , Aug. 1971, pp. 77-89. 

3. Valckx, F. P. G. : Electrostatic Analyzer for the Detection of Fast Neutral Par- 

ticles. NASA TT F- 11458, 1968. 

4. Berkner, Klaus H. ; Pyle, Robert V. ; and Stearns, J. Warren: Cross Sections for 

Electron Capture by 0. 3 to 70 keV Deuterons in Hg, HgO, CO, CH^ and CHF^g 
Gases. Rep. UCRL- 18738, Univ. California Lawrence Radiation Lab. , June 2, 
1969. 

5. Barnett, C. F. ; and Ray, J. A.: A Calibrated Neutral Atom Spectrometer for 

Measuring Plasma Ion Temperatures in the 0. 165- to 10-keV Energy Region. Nucl. 
Fusion, vol. 12, no. 1, Jan. 1972, pp. 65-72. 


10 



60 cm 



Figure 1. - Hot-ion-plasma (HIP-1) source. 


Plasma 



Figure 2. - Neutral particle energy analyzer. 
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Figure 3. - Typical unreduced neutral particle energy spectra 1(E) 
as function of energy. Applied voltage between electrodes, 8 kil- 
ovolts; axial magnetic field at midplane, 1. 35 teslas. 



Energy, E, kV 


Figure 4. - Reduced energy spectra as function of 
energy for applied voltage between electrodes of 
8 kilovolts and axial magnetic field at midplane 
of 1.8 teslas. f(E) ~ I(E)/o 10 E 3 K GC . 
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, arbitrary units 



Figure 5. - Reduced energy spectra as function of AE (atomic number 
times energy) for applied voltage between electrodes of 8 kilovolts and 
axial magnetic field of 1. 8 teslas. Centerline sighting. 



Figure 6. - Reduced energy spectra as function of energy for 
applied voltage between electrodes of p kilovolts and axial 
magnetic field at m idplane of 1. 35 teslas. 



(a) Effect of magnetic field. Applied voltage between electrodes, 8 kilovolts. (b) Effect of applied voltage. Axial magnetic field at midplane 1 8 

teslas. 




(c) IlEjI/Kl keV) as function of magnetic field. (d) IlEjI/Kl keV) as function of applied voltage. 


F 'fo U r r hv 7 drooen rametriC Varia,i ° nS of neutral particle ener ^ analyzer s ‘ )ec,ra I,EI with a PP |i «< 1 voltage and axial magnetic field. Centerline sighting 




Sighting line, cm from centerline 
(a) Hydrogen. (b) Deuterium. 

Figure 8. - Neutral particle energy current 1(E) as function of sighting position. 



A 0 - 486. 1 nm \ 0 - 486. 1 nm A 0 ■ 486. 1 nm 

(a) Sighting 0.9 cent i- (b) Centerline sighting. (c) Sighting 0.9 centi- 
meters above centerline. meters below centerline. 

Figure 9. - Typical H(3 spectra. Applied voltage between electrodes, 8 kilovolts; axial 
magnetic field at midplane, 1. 35 teslas. 



Sighting line, cm from centerline 

Figure 10. - G(A ■ 486. 1 nm) as function of sighting 
position. 
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